Role of beam propagation in Goos-Hanchen and Imbert-Fedorov shifts 
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We derive the polarization-dependent displacements parallel and perpendicular to the plane of incidence, for 

a Gaussian light beam reflected from a planar interface, taking into account the propagation of the beam. 

Using a classical-optics formalism we show that beam propagation may greatly affect both Goos-Hanchen and 

Imbert-Fedorov shifts when the incident beam is focussed. © 2008 Optical Society of America 

OCIS codes: 240.3695, 260.5430. 



It is known that the behavior of bounded beams of 
Ught reflected from and transmitted through a planar 
interface differs from that exhibited by plane waves, the 
latter being ruled by Snell's law and the Fresnel equa- 
tions."'^ For bounded beams diffractive corrections occur, 
the most prominent of which the so-called Goos-Hanchen 
(GH)^ and Imbert-Fedorov (IF) shifts'^ of the beam, oc- 
curring in the directions parallel and perpendicular to 
the plane of incidence, respectively. In principle, both re- 
flected and transmitted beams are subject to such shifts. 
A great deal of literature exists about experimental'*'^ 
and theoretical^"^ demonstrations of both GH and IF 
shifts but generally the effects of beam propagation are 
not accounted for.**^ A notable exception is a recent pa- 
per by Hosten and Kwiat^* where the authors report, 
among other issues, on a dramatic signal enhancement 
technique (^ 100 x) for a quantum version of the IF 
shift, the spin Hall effect of light (SHEL). based on beam 
propagation.*^ The theoretical discussion in** uses the 
quantum formalism of weak measurements,*'^ although 
the authors note that the beam propagation enhancement 
(BPE) is essentially a classical phenomenon. 

The purpose of this Letter is to present a purely clas- 
sical analysis of the BPE; we feel that this is useful since 
a classical description will make this important tech- 
nique, which allows sub-nm sensitivity,** better accessi- 
ble to the metrology community. Furthermore, our clas- 
sical framework covers both the GH and the IF case, 
whereas the treatment in Ref.** is restricted to the IF 
case only. A last, minor, difference between the present 
work and that of Hosten and Kwiat** is that the lat- 
ter authors measure the beam that is transmitted across 
an air-glass interface, while we study the beam that is 
reflected by such interface. Since the transmission and 
reflection cases have a very similar mathematical struc- 
ture, all our main conclusions regarding reflective GH 
and IF shifts remain qualitatively valid for the transmis- 
sion case. 

We begin by considering optical reflection from a pla- 
nar interface; Fig. 1 illustrates the coordinate system. 
The z axis of the laboratory Cartesian frame (xyz) is 



normal to the planar interface (z = 0), that separates 
empty space (in practice air), where z < 0, from an op- 
tically dense region (either a dielectric or a metal*^), 
where z > 0. We use a Cartesian frame {xi,yi,Zi) at- 
tached to the incident beam and another one {Xr, yr, Zr) 
attached to the reflected beam. Note that the coordinate 
Xr is associated with the longitudinal GH shift, while 
yr is associated with the transverse IF shift. Consider 
a monochromatic Gaussian beam of light propagating 
in air parallel to the positive Zi axis. The electric field 
amplitude of such beam can be written as^ 
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where we have introduced the dimensionless variables 
Xi — kxi, Yi = kyi, Zi = kzi, A = kL, and where 
the paraxial approximation corrected up to first order 
derivatives has been used.*^The Gaussian amplitude of 
the beam is characterized by the minimum waist wq lo- 
cated at Zi = 0, and the Rayleigh range L = kwQ/2. 




Fig. 1. (Color online) Geometry of beam reflection at the 
air-medium interface. 



The polarization of the beam is determined 
by the complex-valued unit vector f = 



1 /2 

(/px, + fsft ) / ( I /p P + I /s P) , which corresponds 
experimentally to a polarizer perpendicular to the beam 
central wave vector k^ = fcz^. Across the interface z = 0, 
the tangential components of the electric and magnetic 
field must be continuous. From this boundary condition 
and Eq. (P), the electric and magnetic fields of the 
refiected beam can be determined, under the same 
conditions. A straightforward calculation^ yields: 
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where Xr = kxr, Yr — ki/r, Zr = kzr, and rA,drA/d9i^ 
are the Fresnel reflection coefficients and their flrst 
derivatives evaluated at the "central" angle of incidence 
9i ~ arccos(ki • z/k), respectively. The index A G {P, S} 
is a label for a linearly polarized plane wave whose elec- 
tric field vector is either parallel (P) or perpendicular 
(S) to the plane of incidence {x, z), which is defined as 
the common plane of the central wave vector k^ = kzi , 
and the normal to the interface z. In a similar man- 
ner the magnetic field B''''^ of the reflected beam can be 
obtained and used to calculate the beam intensity spa- 
tial proflle I{Xr,Yr, Zr) as the flux of the time averaged 
Poynting vector S oc Re(E''°^ x B''°^ ) through a surface 
perpendicular to the central direction of propagation z,,: 
I(Xr, Yr, Zr) (X S-Zr- At any given plane Zr — const., the 
intensity I{Xr, Yr, Zr) can be considered as the distribu- 
tion of the "quasi-Gaussian" variables X = XrXr + YrYr 
with means M = {Xr) Xr + (Yr) Yr, where 
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determines the centroid of the reflected beam.^^ As the 
result of a straightforward calculation, one obtains for 
the GH shift 
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and for the IF shift 
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where r^ = i?ACxp(i0A), A E {P,S}, fp ~ ap E M, 
fs = a5cxp(i?7), and pA = Rc{dlnrA/ddi), (pA = 
Im {d\nrA/d6i). 

Equations (|4l5p are the first main result of this Letter. 
They give both the GH and the IF shift as functions of 
the beam propagation distance Zr- For a well coUimated 
beam, the condition Z^/A ^ 1 is trivially satisfied at op- 
tical frequencies and the Zr-independent terms in both 
Eqs.(j4][5|) are dominant. Such terms represent the "tra- 
ditional" GH and IF shifts. In fact, when Zr = 0, Eq. 
([4| is a straightforward generalization of the well known 
Artmann formula; ^^ and Eq. ([5]) is in agreement with 
Bliokh and Bliokh.^^ However, for a focussed beam the 
condition Zr/A ^ 1 may hold, and the Z^-dependent 
terms become relevant. For example, if a typical He-Ne 
laser operating at wavelength A of 633 nm with a mini- 
mum waist wi of 1.5 mm, is focussed by a lens with focal 
length / of 100 mm, a waist wq = Xf /{irwi) ~ 13 pui is 
produced. If this beam propagates over a distance of 250 
mm from the lens, we easily obtain Zr/A ~ 168. Thus, 
depending on the actual experimental conditions, either 
Zr-dependent or Zr-independent terms in Eqs. (|4|5p may 
be dominant, and, as a consequence, the measured GH 
and IF shifts will dramatically change. 

In order to show the connections between the results 
above and the ones presented in Ref.,^^ we must take 
a step backward and calculate, as an illustrative exam- 
ple, the quantity x^ • E''°^ evaluated for ag = on the 
plane of incidence Yr = 0. This corresponds to an exper- 
imental configuration apt to measure the GH shift of a 
P-polarized beam. From Eq. ([2]) it readily follows 
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A careful inspection of the equation above, shows that it 
represents a Gaussian beam tilted clockwise by an angle 
Pp/A with respect to the axis z^, and displaced by ipp/k 
along the axis x^. In other words, for a Gaussian beam an 
imaginary position shift is equivalent to an angular tilt, 
while a i^eal one represents a spatial shift. It is easy to 
see that the complex- valued nature of the shift, controls 
its behavior under beam propagation. To demonstrate 
this, we note that calculation of the centroid of the dis- 
tribution |x,. • E^^p yields 



(Xr) =PP- {Zr/A)pi 
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which coincides with Eq. ([4]) evaluated for a^ = 0. This 
expression of {Xj.) shows that only the imaginary part pp 
of the complex shift Lpp — ipp couples to the coordinate 
Zr and it is enhanced by a factor Z|./A as the reflected 
beam propagates along the optical distance Z^. As it 
will be shown below, such factor Z^j h. coincides with 
the propagation enhancement factor F of Ref.^^ which 
lies at the core of their signal enhancement technique. 
Although in the reasoning above we have considered a 
specific example, it is not difficult to realize that the 
conclusions reached are perfectly general; in particular, 
positioning of the waist of the Gaussian beam at the 
interface is not essential (see also Ref.^^). This is the 
second main result of our Letter. 

We conclude this Letter by demonstrating that, as an- 
ticipated, the term Z^/A in our Eq. ([7]) is coincides with 
the propagation enhancement factor F of Hostcn and 
Kwiat.^^ According to their scheme, the incident beam 
is first pre-selected in the P polarization state, namely 
/p = 1, /s = 0, and then post-selected (after reflec- 
tion by an air-glass interface) in the v polarization state, 
where v = x^sinA -I- y^cosA, and A = ±|A|. In ad- 
diction, after reflection the beam is observed along the 
transverse plane X^. = 0. Thus, the relevant amplitude 
of the reflected field can be written as 
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where the purely imaginary term iD cot A is described in 
Ref.^^ as the product of the SHEL-induced photon dis- 
placement D = {1 + rs/rp) cot 9i, and the weak value 
of the photon spin component i cot A. However, by com- 
paring Eq. (HI) with Eq. (I6|) , it is clear that at a classical 
level such imaginary shift amounts to a tilt by the small 
angle D cot A/A, of the reflected beam. From the dis- 
cussion above we know that an imaginary shift couples 
with coordinate Zr and, therefore, increases as the beam 
propagates. In fact, calculation of the centroid of the dis- 
tribution |v.E''''^|^ ^^ p yields 

{Yr.) = iZr/A)DcotA, (y,2) = (a2 + Z2)/(2A), (9) 

where Zr/A = i^ (~ 156 in Ref.^^) is just the propaga- 
tion enhancement factor. This last equality can be easily 
proved by combining the two formulas in Eq. ^ to ob- 
tain, in terms of dimensional variables, 
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where S = D/k, R{z) = {z'^ + L'^)/z 



is the ra- 



dius of curvature of the beam wavefront at the lens L2 
position, ^^ and the last approximate equality holds in 
the experimental conditions of Ref.,^^ where z/L ^ 1. 
From Eq. ^ it readily follows that F = 47:{yl)/{z\), 
in agreement with Eq. (4) in Ref.^^ 



In summary, we have furnished analytic expressions, 
based upon classical optics, for both the GH and the IF 
shifts, as functions of the beam propagation distance. 
These give in a natural way the dramatic changes of 
the GH and IF displacements induced by using a fo- 
cussed beam. Moreover, from the analysis of such ex- 
pressions, we derived a fully classical interpretation of a 
very recently introduced signal enhancement technique, 
employed to measure the spin Hall effect of light. ^^ 

We acknowledge Michele Merano for useful discus- 
sions. This project is supported by FOM. 
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